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We present here how to study steady hnear transport phenomena by using entropy. 
We study the system and the environment together and identify their entropies. 
Concerning their interaction, quantum mechanics is considered in the study. A 
time parameter r is therefore introduced to characterize the discrete nature of the 
quantum interactions. By combining r and the entropy, an approach is constructed 
successfully to study electrical conduction, heat conduction and shear viscosity. 

PACS numbers: 05.60.Cd, 05.70.Ln, 05.70.Ce, 44.10.+i, 02.50.Cw 

Keywords: transport, entropy, heat conduction, electrical conduction, viscosity 

INTRODUCTION 

In this article, we study three hnear transport phenomena: heat conduction, electrical 
conduction and shear viscosity. All the three phenomena are linear when the external applied 
forces are weak. These linear behaviours are described by empirical laws. Heat conduction 
is described by Fourier's law 

3 = -kVT, (1) 

where j is heat current density, T is temperature, k is thermal conductivity. Electrical 
conduction is described by Ohm's law 

j = aE = -aVU, (2) 

where a is electrical conductivity. Shear viscosity is described by Newton's law of viscosity 

= (3) 

where Tzx is shear stress, rj is shear viscosity. 

The three laws are in similar forms. So they are likely associated with entropy which 
is the most general concept. About a system in the equilibrium state, its entropy is well 
understood. But about a system in a non-equilibrium state, its entropy is even not well 
defined. This is because the former is static while the latter is dynamic. The three transport 
phenomena belong to non-equilibrium state and therefore their empirical laws have not been 
theoretically derived. Yet, the three phenomena are steady and linear. Their dynamical 
behaviour should follow a very simple principle. In this article we will use entropy to 
construct this principle. 

The present article is organized as follows. In section 2 we review the entropy of the 
equilibrium state. In section 3, we try to apply the entropy to steady non-equilibrium state 
and identify the problem. In section 4, we solve the problem by quantum mechanics. In 
section 5, we study electrical conduction. In section 6, we study heat conduction. In section 
7, we study shear viscosity. 
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OVERVIEW OF ENTROPY OF EQUILIBRIUM STATE 



We consider the system and its environment as two parts. Let's denote the total energy 
by Eq and energy of the system by E. The environment energy is then Eq — E. For a given 
E, the system (environment) can be in a number of possible states. That number is denoted 
by Q {Q') from which the entropy is defined as S = klnQ {S' = k\nQ'). Thus we have 
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The environment should be very large and therefore its entropy linearly depends on the 
energy fluctuations. So the temperature can be introduced by 
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(6) 
(7) 

(8) 



which is actually the Boltzmann distribution. The average of E can be calculated through 
the partition function 

Z = J2^iE)e-^ (9) 

E 

which leads to 

1 

(10) 



E = ^^Efi(E)e~^. 



The most probable value of E can be found by maximizing Eq. 
logarithm 

which is equivalent to minimizing the free energy 

E-TS{E). 



or maximizing its 
(11) 

(12) 



The most probable E and the average E will be equal and will be actually the observed 
energy when the system is big. But it should be noted that when the system is small, they 
are not necessarily equal and also the observed energy will fluctuate explicitly [sj . The same 
argument applies to for example electric current. 
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STEADY NON-EQUILIBRIUM STATE 



We study here the steady non-equihbrium state by an example of an electrical conductor. 
Let an electric field turn on at time t = 0. The electric potential is U. Let the electric 
current I be steady. Then at time t, it will generate energy Ult which is absorbed by the 
conductor and the surroundings. We can simple take the whole conductor as a system. But 
we had better take only those movable electrons/ions as a system. Then the surroundings 
and the rest of the conductor such as the crystal lattice together form the environment. 
Anyway, at a time t Eq. @ becomes 
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(13) 



and Eq. ([8]) becomes 



P{EJ) oc J)exp 



Ult-E^ 



(14) 



When the conductor is big enough, the most probable current and the average of the current 
are approximately equal and are actually the observed steady current. We will only calculate 
the most probable current because the calculation is easier. The most probable current can 
be obtained by maximizing P{E,r). To do that, we need to separate / from E. So we 
assume 

Q{E,I) = Q{E)Q{I). (15) 

Its meaning can be understood by an example. Consider a crystal lattice and two interstitial 
ions. The ions can jump from site to site beside vibrating. Suppose the vibrating energy is 
one of {0, hcu, 2?ilo, ■ ■ ■}. Also suppose each ion jump once and each jump produces electric 
current proportional to either qa/r or —qa/r. Then for E = 2%uj and / = 0, we have 



fi(2^w,0) = 3x2. 



(16) 



Thus Eq. f ll4p can be factorized into two parts. One part only contains E, the other part 
only contains /. The E part will lead us to the free energy and subsequently the most 
probable E. The / part we hope can lead us to the most probable current. This part is 



P(/)ocfi(/)exp(^). 



Let's introduce system entropy 



S{I) = k\nn{I). 

Then it seems that the most probable / can be obtained by maximizing 

Ult 



Sil) + 



T 



(17) 
(18) 

(19) 



But this is impossible because t is a variable. The problem arises because we have overlooked 
the quantum nature of the microscopic interactions. 
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QUANTUM MECHANICS AND r 

According to quantum mechanics, a system evolves unitarily. But the unitary evolution 
is constantly interrupted by the interactions coming from the environment. The interaction 
induces the system wavefunction to collapse from time to time. Suppose this happens every 
r. Then during time < t < r, the system is in a deterministic wavefunction and the 
wavefunction evolves unitarily. But at time t = t, the wavefunction collapses and the result 
of the collapse follows a probability distribution from which the entropy emerges. Thus Eq. 
(fT3|) only makes sense at time t = t and other moments that the wavefunction collapses. So 
if we redo the calculation about Eq. (fT9l) . we will get 

UIt 

S{I) + -jr. (20) 

from which the most probable / can be solely determined. To understand it better, we apply 
it to an example. 

Let's consider here the ionic conduction in which case the conductor is a crystalline solid 
conductor that is composed by a non-movable crystal lattice and some movable interstitial 
ions. An ion can jump from site to site. Suppose each ion jumps once every r and each 
jump produces electric current proportional to ±ga/r, where q is the ion charge, a is the 
lattice constant. Suppose the ion concentration n is very low. 

First, let's consider a system that is composed by just one ion. And let the electric field 
E turn on at time t = 0. Then at time t = t, the ion will either jump up t or jump down 

So the result will be 



state : { | t) | |) }, 

ly = { qa/r -qa/r }, 

Q = { 1 1 } 

P oc { 1 X exp(ff ) 1 X exp(-|f ) }, 



(21) 



where 



Iy = J2qv,. (22) 

i 

But during time < t < r, the system is a pure wavefunction 

\^^I)=A^\t)+A2\i), (23) 

where l^ip cx exp{j^), |y42p oc exp(— ^). Only at time t = r, it will collapse. When it 
is going to collapse, it can somehow sense the number of the final states and tune |y4i| and 
1^421 to the right amplitudes. The number of the final states for this example is the number 
of the environment states. But that number depends on the result of the collapse that does 
not happen yet. If the result is | ']'), the environment will gain an extra energy Eaq and 
thus gain extra entropy Therefore that number will increase by a factor of exp(-||?). If 

the result is | J,), the environment will loss energy Eaq and thus loss entropy Therefore 
that number will decrease by a factor of exp(— ^|^). Once the wavefunction has collapsed, 
another unitary evolution will start and the result will be Eq. f l23|) again just before next 
collapse. About how a wavefunction collapses, it is a quantum measurement problem which 
is studied by the decoherence theory @, [7|. 
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For a system that is composed by two ions, Eq. (12T]) becomes 

state: { | t)| t) | t)l i) or | i)| t) | i)| i) }, 

Iv = { 2qa/T -2ga/r }, , . 

n = { 1 2 1 }, ^ ' 

P oc { 1 X exp(^) 2x1 1 X exp(-^) }. 

For a system that is composed by ions, it becomes 

Iv = { Nqa/T ■■■ {2k-N)qa/T ■■■ -Nqa/r }, 

n = { 1 ■■■ ■■■ 1 }, (25) 

P oc { 2xexp(^) ■■■ X exp((^^^^) ■■■ 1 x exp(-^) }. 

When N is large, we can use an approximation of 

C^-2-^expf-i^l^). (26) 

Then by writing k into ly and by I = ly/h, we get the system entropy 

- ^» - il^^ ■ '''' 

where 5*0 is the entropy of the equihbrium state, h is the height of the conductor. Thus Eq. 
( 12 up in this example becomes 

S(I) + — = S. - -^/^ + — , (28) 



by which the most probable / is obtain as 



krhT 

By the way, the electrical conductivity is also obtained 



E. (29) 



which is the same as the result of the kinetic theory 

One can also calculate the most probable / directly from Eq. ( |25ll . But our calculation 
by Eq. (I2U]) is easier. Besides, Eq. (12U]) is in a general form and is ready to be applied to 
other studies. 

We see that r is actually the time interval between two successive jumps of the same 
ion. About the time interval between two collapses of the whole system, it is much smaller 
especially when the system is large. For a system containing N ions, its wavefunction 
collapses approximately times during r. We also see that r is actually an average quantity. 
If we at the beginning let ions jump randomly, we will still obtain all the results the same. 
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ELECTRICAL CONDUCTION 

We have studied a specific electrical conduction. Now we study the general electrical 
conduction. Thus we need to write the system entropy in a general form like 

Sil) = So- CJ - ^C2l' - ^Cil' + ■■■. (31) 

5*0 is the entropy of the equilibrium state. Ci should be zero because the point of / = 
should be a stationary point. Besides, S{I) should be an even function. Therefore when / 
is small, the system entropy can be approximately written as 

3(1) = So- \C2P. (32) 

We have not specify the shape of the conductor. C2 actually depends on it. Let's introduce 
the system entropy per unit cube as 

s{i) = So- ^C2j^ (33) 

where j is the current density. Then for a conductor having cross section A and height h, 
there will he S = sV, I = jA, C2 = C2h/A. 

For the ionic conduction, we have C2 = ^^2^2 which is proportional to r^. So C2 is actually 
not independent. We had better introduce an independent parameter 

b = ^. (34) 

By b, we can recast s{j) and we will see that it is actually independent of r. It is actually a 
configurational entropy which is about what the probability for a certain kind of microscopic 
configuration to appear. 

)]) becomes 

-^C2j (35) 

j = 7^E. (36) 
brT ^ ^ 

This formula is the same as Ohm's law ([2]) . In addition, the electrical conductivity is 
obtained. 



Concerning j, Eq. (l20l) becomes 



which leads to 



HEAT CONDUCTION 

The heat conduction can be studied similarly. Let's take a system and fix it in between 
two thermal baths. One bath has temperature Ti, the other T2. Let Ti < T2 and let a 
steady heat current J flow through the system. As an analog of Eq. (133|) . we have 

sU) = So- \c23\ (38) 
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And as an analog of Eq. (132|) , we have 

1 



S{J) = So- ^C^J^ (39) 



Thus about heat conduction, Eq. (l20l) becomes 



lev' + Jt{^ - (40) 



which leads to 

Let the system have height h. In the limit of /i — )■ 0, the above equation becomes 

J--^^ (42) 

^~ C^T^dh ^^^^ 

which can be recast as 

where h = c^jr^ which is introduced as an analog of Eq. This is the same as Fourier's 
law ([1]). In addition, the thermal conductivity is obtained. 



For dilute gas, has been calculated in a paper [9|. By using Eq. (HT]) . the paper has 
also calculated the thermal conductivity and the result is the same as the kinetic theory. In 
that paper, Eq. (HT!) is obtained by an analysis of the entropy production. 



SHEAR VISCOSITY 

The previous analysis can be applied to shear viscosity too. Think two 2;-axis perpendic- 
ular plates between which it is filled with a fluid with cross section A and height h. Let the 
bottom plate keep still, the top plate move along x-axis at a steady speed Vx- Each molecule 
of the fluid takes part in two kinds of movements. One is the heat movement. The other 
is the macroscopic movement associated with The former contributes entropy to the 
environment. The latter contributes entropy to the system. 

Let the velocity gradient ^ be the same everywhere. Then as an analog of Eq. (!33l) . 
the system entropy per unit cube is 

Avx. 1 ,dvx^2 

^< d7' = - ■ 

During r, the force F applied to the top plate will generate heat Fv^t which subsequently 
brings an extra entropy Fv^t/T. Then as an analog of Eq. fl35|) . the steady velocity gradient 
can be obtained by maximizing 

1 ,dVx.o TzxT dVx s 
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which leads to 

where Tzx = Fj A. It can be transformed to the form of 

r..^..T^. (48) 

where h = C2/t'^ which is introduced as an analog of Eq. flM|l . It is the same as Newton's 
law of viscosity ([3]). In addition, the shear viscosity is obtained, 

r] = brT. (49) 



CONCLUSION 



We have studied electrical conduction, thermal conduction and shear viscosity by the 
same approach. In the approach, we have used the entropy as a central concept. We have 
also considered quantum mechanics explicitly and introduced a parameter r to characterize 
the time interval between two successive interactions of a microscopic particle. This allows 
us to analyze the system/environment interactions as a series of discrete quantum processes 
from which entropy emerges to characterize the probability for a microscopic configuration 
to appear. We have taken the steady electric current as the most probable current and 
established an equation to find it. The found current follows Ohm's law. In the same way, 
we have reproduced Fourier's law and Newton's law of viscosity. 
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